Comment on 'Ising Pyrochlore Magnets: Low Temperature Properties, "Ice Rules," 

and Beyond' by R. Siddharthan et al. 
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Siddharthan et al. discuss the competition be- 
tween dipolar coupling and superexchange in pyrochlore 
magnets with (111) Ising-like anisotropy, such as 
Ho 2 Ti 2 7 |] and Dy 2 Ti 2 7 §|. In the simplest ap- 
proximation for Ho 2 Ti 2 07, with Ising spins and nearest- 
neighbor ferromagnetic exchange, one obtains the "spin 
ice" model that predicts a macroscopically degenerate 
ground state S. Consistent with this, /xSR and neutron 
scattering data find no phase transition in Ho 2 Ti 2 07 
in zero field down to at least 50 mK, and at that temper- 
ature the neutron scattering pattern is consistent with 
the spin ice ground state p). On the other hand, the 
dipolar coupling, Jd, in these materials, is of the same 
order of magnitude as superexchange, Js- Siddharthan 
et al. derive an estimate for J5 = — 1.92 K for the 
near neighbour exchange and Jd = 2.35 K. These val- 
ues are consistent with our own estimates Jd — 2.35K, 
and Js = —1.4 ±0.5 K from single crystal measurements 
where demagnetizing effects were explicitly accounted 
for, and which correspond to the Curie- Weiss tempera- 
ture of +1.9 K, quoted in Ref. Using their estimated 
value of Js, Siddharthan et al. simulated a model in 
which the dipolar summation is truncated at five nearest 
neighbours, and find a transition to a partially ordered 
phase marked by a sharp peak in the heat capacity, C (T) , 
at T ~ 0.8 K Q. They further argue that such a pic- 
ture applies to Ho 2 Ti 2 07 by commenting that the exper- 
imental heat capacity shows, as their simulation does, a 
steep increase at about 1 K and that the experimental 
data for Ho 2 Ti 2 07 suggest a vanishing entropy for the 
ground state (on integrating C(T)/T). Their conclusion 
of a phase transition in Ho 2 Ti 2 07 is inconsistent with 
the experimental /LtSR and neutron data of Refs P,H,@] • 

In view of the very large hyperfine coupling in Ho, as 
observed in many other Ho-containing compounds such 
as H0F3 H , it is necessary to include a component from 
the nuclear spin, /. Blote et al. M find a heat capacity 
in the pyrochlore Ho 2 GaSb07 very similar to the one 
of Ho 2 Ti 2 07, which they model below 2 K as a simple 
Schottky anomaly with the theoretical maximum for Ho 
(/ = 7/2) of 0.9 R at a temperature of about 0.3 K. 

In this Comment we argue that Ho 2 Ti 2 07 does in fact 
exhibit spin ice behavior, and that the experimental spe- 
cific heat data can be accounted for in terms of a dipolar 
spin ice model || and by including an expected contibu- 
tion from the nuclear spins |Q to the appropriate long- 
range treatment of the dipole-dipole interactions. 

Dipole-dipole interactions are conditionally convergent 
due to their 1/r 3 nature and their lattice summation 
must be considered with care Q|. In order to include 
the long range nature of the dipole-dipole interaction, 
we have used the standard Ewald method in our Monte 
Carlo simulations with either of the above set of {Js, Jd} 
values. We found it sufficient to simulate 4x4x4 cubic 


cells (1024 spins) with ~ 10 6 Monte Carlo steps per spin. 

Figure 1 shows the electronic spin part heat capacity 
for J s — — 1 K and Jd = 2.35 K (open squares). The to- 
tal entropy is within 2% of i?{ln(2) - (1/2) ln(3/2)} and 
therefore agrees with that of spin ice . We have found 
in our Ewald simulations that dipolar spin ice occurs for 
all values Js/Jd > —0.91, which includes the value of 
Rcf . jl] , and that fully developped Q — long-range an- 
tiferromagnetic order occurs for Js/Jd < —0.91, with 
no partially frozen or ordered state as found in Ref. [|lj. 
The open circles show the sum of this magnetic heat 
capacity with a nuclear Schottky anomaly as found in 
Ho 2 GaSbC>7 @ . We find that these results reproduce rea- 
sonably well the experimental data presented in Ref. |l[] 
(filled circles), and is definitely more accurate than the 
model proposed by Siddharthan et al. Equivalently, sub- 
tracting the nuclear contribution from the magnetic heat 
capacity measured in Ref. Q leads to a heat capacity sim- 
ilar to that found in the spin ice material Dy 2 Ti 2 07 [[|. 
In conclusion, we argue that the experimental specific 
heat of Ho 2 Ti 2 07 can be accounted for by the sum of a 
dipolar spin ice magnetic contribution and a nuclear hy- 
perfine Shottky anomaly with no need to invoke a tran- 
sition to a partially ordered state. 
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FIG. 1. Open squares: heat capacity of dipolar spin ice 
with J s — — 1 K and Jd = 2.35 K. Open circles: same data 
with the additional nuclear contribution expected for Ho [6,7]. 
Stars: experimental data extracted from Ref. [1]. 
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